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I. GENERAL FORMALISM 



The quark parton model of Feynman [l| provides the simple description of deep inelastic phe- 
nomena as well as Drell-Yan processes. Theoretical justification of this model was done in terms 
of asymptotically free gauge theories 0| . Drell-Yan picture based on factorization of contributions 
from small and large distances was justified in papers of Collins Deviation from the naive 
Bjorken scaling of the structure functions of deep inelastic scattering (DIS) was recognized to be 
broken by the so called "large logarithms" — the logarithms of the ratio of momentum squared 



Q 



-q (virtualities) of particles deep off-mass-shell to their masses. The reasons of appearing of 



these logarithms in QED was clarified by methods of quasi-real photons and electrons [j-Ql ■ Keep- 
ing in mind the contributions of higher orders of perturbation theory,in the leading logarithmical 
approximation description of processes with large virtualities can as well be formulated in parton 
language with some definite dependence on of partons (quark, gluons) densities — structure 
functions, 9' {x,t), G{x,t). 

The evolution equations of Altarelli-Parisi (AP) [8J for these densities have a form 



dq^ [x, t) _ as jt) i dy 
dt ~ 2Tr J y 



G{y,t)P,G 



t = In 



(1) 



dC ix,t) ^ as jt) rdy 
dt 2n J y 



2Nf 

J2q'{y,t)Pa,.('^j+Giy,t)PGG i^^ 



, i ~ u,d,s, c, b, t. 



They describe the dependence of quarks {q{x,t)) and gluons {G{x,t)) densities with the energy 
fraction x inside the proton on a scale of (where > 0). Here m is the suitable normalization 
point m ^ Qq ^ Mp 1 GeV, as (t) is the QCD coupling constant on the scale and Pij {x) 
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Figure 1: Definition of correlator densities. 



are the AP equation kernels: 

with Cf = and = 27V for the color group SU{N). These quantities satisfy the following 

properties: 



dz P 



i 

(z) =0, Jdzz Pgg (z) = 0. (2) 



It's useful to remind here the statistical interpretation of AP equations in terms of densities @ 
by means of a set of correlation functions, satisfying the system of statistical equations (renormal- 
ization group equation). It was a success of the numerous applications of APL set of equations, 
working with two densities q, G of quarks and gluons into a proton. 

Problems associated with processes with large multiplicity [lol - [l2j , however require some gener- 
alization of traditional approach, introducing the correlation functions . 

Namely let us introduce D'^{x,t), D3{x,t), D^{x,t) as the densities of quark, gluon, anti-quark 
into the initial quark and the similar quantities for the initial anti-quark (see Fig. [1]). 

Besides let us introduce G''{x,t), G^{x,t), G^{x,t) as the similar densities for the initial gluon. 

In complete analogy to the case of APL equations one can obtain the evolution equations for 
these set of densities. They are presented in Appendix \K\ 

When neglecting the presence of densities of anti-quarks in gluon and quark, the combinations 
G^{x,t) + D^{x,t) = G{x,t) and G^{x,t) + D^{x,t) = q{x,t) can be shown (see Appendix E| to 
obey the equations of AP. 

One of successful phenomenological model-based on the dominant role of gluon distribution 
in describing the processes with high hadron multiplicities [10 12]. 

Solving the equation for G^ by iteration method we see that regeneration of gluon density in the 
channel G^ — 13' — > G^ is associated with small factor 

Terms of such a magnitude can be neglected thus determining the accuracy of the approximation. 
Alternatively it can be included as a some contribution to K -factor. 

Let now identify the nonsinglet structure function q — q with D'' - a quark dominance density. 
We introduce the gluon dominance density G^, satisfying the equations: 
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s 1 

xG^x,/?,) = Six-1) + J J ^PGG{^^)yD,{y,P'^), (5) 

where s is the total invariant mass of the process (i.e. ^/s = 2E, where E is the energy of particle 
in beam) and 



2 



Using the solutions of the homogeneous equations for quark non-singlet density D'^ and G^: (we 
use the method similar to one developed in frames QED in (see Eq. (11) and Eq. (20) in [14]): 



xG^ {x, /?,) = 2/3g (1 - xf^^ (^1 + ^/3,^ -2x + x^ (1 - x)^ + O (/J^) 



(7) 



These solutions satisfy the properties i.e.: 

1 1 

dx D'^ {x,l3g) = 1, I dx xG^x,l3g) = 1. (8) 







Below we consider the Drell-Yan process in collision of protons with the hard subprocess a+b — >■ 
(where a and b are partons and F"^ is some final state produced by them) which is the part of 
more complicated process p + p ^ jeti + iet2 + F"-^ . Thus we choose scale of order = s (i.e. 
O-s (Q^) = cts) where y/s — 2E is the total energy of process in the center of mass frame. 

For inclusive experiments in pp-coUisions the cross section in center-of-mass system takes form 
of Drell-Yan type form: 

1 1 

dapp-^F+x = j dxi J dyi^Wai{xi)^Dl\{xiyi,Pb,) x 
1 1 

X / y dy2 Walx2) ^ Dl\{x2y2, Pb^ ) X 

X d&''^''^^''{s,i,u)Q{z- zth) , z = xiyiX2y2, (9) 

where W"" {x) is the probability to find a parton of sort a with energy fraction x inside a proton 
with small virtuality (module of it's momentum square of order of 1 GeV^), these quantities was 
obtained in [l^ (see Appendix [Ct . D'^ {x,f3) is the densities of parton of sort b inside the parton 
of sort a. The summation over {01,02} is performed over all possible partons inside the proton, 
i.e. It, d, s, u, d, s, g. The summation over {61,62} is performed over all possible partons which 
can be found inside the parton of sort {ai, 02}, i.e. in principle all possible partons too (u, d, s, u, 
d, s, g). The quantity 

da''^b^^^{s,i,u) (10) 

in is the cross section of hard subprocess of two partons 61 and 62 fusion which actually produces 
the final system F which is of experimental interest. This cross section is considered already in 
the system of center-of-mass of this subprocess 61 4- 62, i.e. these invariants |s,t,M} are in this 
reference frame of subprocess and have a form: 

s = AE'^z, i= -2E^z(i-cos9) , u = -s-i, 
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where z = X\y\Xiyi and d is the angle between 3-vectors of initial parton h\ momentum and the 
momentum of one of particles from the created state F in the center-of-mass reference frame of 
the subprocess, which can be expressed in terms of the angle Q between the directions of the initial 
beam and the momentum of one of particles from the created state F: 

cos6' = -. (11) 

X2y2 + xxy\ cosO 

Here polar angle 9c is the angle between 3-vectors of initial parton A and the momentum of one of 
particles from the created state Fab in center of mass frame of sub-process. Polar angle 9 is the 
angle between the directions of the initial beam and the momentum of one of particles from the 
created state Fab- 

The quantities Ka in ^ are the so called i^-factors which takes into account non-leading con- 
tribution of evolution. The if-factor associated with quark density has a form p^ : 

1 , , , 1/67 TT^ \ 5n/ 
if o = 1 H fco, ka = -Cv ~ 1.5. 

and Kg is the iiT-factor associated with gluon density and has a form 

1 
2 



Kg^l + -Ko (12) 



where Kq is given in ([3]). 

The 0-function in ([9]) assures that experimental setup allows to register the jets of produced 
particles with some finite invariant mass Sth only, i.e. the jets with the invariant mass Sj > Sth- 
And the quantity zth which characterizes this threshold has the form: 

zt,. = — . (13) 
s 

The solutions of evolution equations for = D'^, are given in ([7]). The solutions for = 
, D^, — are presented in Appendix [B|) . Density correlators of type 13' which describe the 
density of a quark q' in quark q and (the density of anti-quark q' inside the quark q) as well 
are discussed in Appendix [Bl 



II. APPLICATION TO SOME DEFINITE SUBPROCESS 

Below we consider two types of subprocesses. First, for the sake of demonstration we will consider 
the process of of associative production of top quark and the charged Higgs boson since the 
application of our approach is more simple in this case. And then we will use our approach to 
describe the experimental data from Tevatron (22j . 



A. The process p + p ^ t -\- H + jijj 

Let us consider now the important application of our approach to the process 

P + P^t + H- (14) 

where j denotes jet. In this case the dominant channel of charge Higgs production is through the 
subprocess: 

b + g^t + H-. (15) 
The cross section of this sub-process has the form (see (2.1) in [3): 

da^g^iH- ^ r J + { _ Mjj_ ml {u - ) + (£ - m^) + s {u - m^) 



2s s (ii — TO( ) 

(u - - s/2) +su/2) 

{u-mlf J' 



(16) 
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where the sub-process invariants s, i, u are defined as: 



S={Pb+Pgy: i={Pb+Pt)\ U= {ph+PH- f , (17) 

and angle 9 is the angle between the momenta of initial 5-quark and produced i-quark in the 
reference frame of center of mass of subprocess (i.e. pb = — Pg where pb and Pg are the 3-momenta 
of initial 5-quark and the gluon correspondingly) ; the quantity do is the following constant 

naos [ml tan^ /3 + cot^ /3) 
~ GM^sin^ 9w ' 

where 9w is the Weinberg angle and /3 is the parameter of Minimal Supersymmetric Standard 
Model (MSSM). For tan^ = 40 and = 0.1 we have ctq ~ 0.06. 

The application of our master-formula ^ to the process gives the cross section in the 
following form: 

1111 



f f f f da''s^*^ 
^ dxi dyi / dx2 / 9 (z - zth) ^^^^q ^ 



dcos 9 

X {Wu {xi)D'i\xiyi,liq) Kq + Wd (xi) Z?«'(xiyi, /?,) + Wg (xi) G«(xiyi,/3g) _ 
X ( Wu (x2) D3{x2V2,liq) + Wd (X2) DS{x2y2,l3<j) + Wg (2:2) G» (2:21/2, /3ff) Kg^ . (19) 

At this stage we need to notice that the cross section of sub-process in (IT51) depends on angle 9 
between the direction of momentum of produced i-quark and the beam in the system of center-of- 
mass of initial proton-proton beams, while the expression (fT6)) depends on the scattering angle 9 
in the center-of-mass reference frame of sub-process. Since these angles correspond to each other 
with the relation (fTTj). i.e.: 

dcos^ = ^yl^AyL_d cos 9, (20) 

{x2y2 + xiyi cos 6*) 

we obtain the following final form of the cross section ([T9| : 

d a„^,H-,,,,, ]dy2 9{z-z y'^'" Avl - Avl 



dcos9 J J J J dcos9 (x2y2 + xiyicos9) 

0000 

X (Wu (xi) D^'{xiyi,l3g) Kg + Wd (xi) 7^«'(xiyi, Kg + Wg (xi) G«(xiyi, /?,,) Kg^ x 

X ( Wu (x2) D3ix2y2,(3q) Kg + Wd (^2) D^X2y2,/3q) Kg + Wg (^2) G^ix2y2, /3g) Kg'^ . (21) 

The dependence of this cross section on scattering angle 9 is presented in Fig. [2] via quantity 

p^^0)^J_^^pp^lEl±2221^ (22) 
do d cos 9 

which is built for few values of charges Higgs boson mass M^. The total cross section of this 
process is proportional to the quantity Th- 



Th{s)= I d9FH{9), (23) 

which is presented in Fig. [3l 
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Figure 2: The angular dependence of quantity Fh defined in (|22|) 
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Figure 3: The quantity Th defined in (|23[l as function of total invariant mass s. 

B. Six-jets production at Tevatron 

In paper the data for six-jet production are present. The application of master-formula ^ 
to this process gives more complicated result since we need take into account few subprocesses [17| : 



d(T{qq qq) _ 



dcosi 



9s 



da{qq' qq') t + w 

dcr,?,0 9s s2 
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Figure 4: The angular dependence of quantity Fqq defined in (|25|) for different initial total energy ^/s. 



dcr{gg qq) _ a. 



dcos ( 



24s ^ 



1 

tu 
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(24) 



where 6 is the angle between the direction of motion of initial parton and the momentum of final 
quark and invariant s, t, u are defined in the same manner as in ([T7|) . 

The comparison of angular distribution of jets momenta in proton-proton scattering from (1241) 



1 da> 



(Ttot d COS 9 



(25) 



(where atot is the total cross section of 6 jet production) with the experimental results (see figure 
6 (a) in paper is shown in Fig. |4] for different initial total energy ^/s. 



III. CONCLUSION 



In this paper we discuss some modification of the method of taking into account the QCD 
leading logarithm radiative corrections based on the Structure Functions approach. Modification 
consists in construction of set of evolution equations for density of parton of sort a in the initial 
quark D° {x,l3g) and density of parton of sort b to be in the initial gluon {x,l3g). This set of 
equations is solved in quark and gluon dominance approximation D"^ :s> D°- , a ^ q and » G", 
a ^ g. This approximation can be improved for the accuracy level which is required by using the 
iteration procedure. This assumption is known as a gluon dominance which is used in description 
of multiplicity of 7r-mesons in the hadron collisions [l(]| - [l2| . 

We present approximate solution for and demonstrate the application of this function 

to the problem of calculation of QCD radiative correction calculation in some particular processes. 

Some efforts was paid to the problem of calculation of subprocesses cross section in the next -to- 
leading approximation. Main attention was paid to the 2-loops level contributions in [l^, [2^ . As a 
result the terms of order (a^L^), (a^L^)^ was taken into account. However the emission of real (soft 
and hard) gluons with the 1-loop radiative corrections was not considered. The role of radiative 
(virtual and real) corrections leads to the change of the (aL^) "-regime to a single-logarithmical 
regime (i.e. only terms ~ {otL)^ remains) in inclusive experimental approach. Single-logarithmical 
approach is determined by renormalization group evolution equation and thus allows us to use the 
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Structure Function approach to obtain the cross section in leading (i.e. (aL)") and next-to-leading 
(i.e. a (aL) ) approximation. 

Functions W°'{x) describe the probability to find parton a inside a proton with off mass shell 
about one GeV squared. These functions was builded in as a result of self-consistent analysis 
of many sub-processes, and besides shown to satisfy the momentum and number sum rules. For 
experimental setups with the product of subprocess detected which moves at large angles with 
invariant mass square exceeding some threshold value sth = s zth the role of " sea" -partons in the 
proton can be neglected. 
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Appendix A: Full set of evolution equations 



Let introduce three distributions Z?" = D^{y,t) with a — q,g,q which describe the number 
of partons a with energy fraction y inside the parent quark. Similarly one must introduce three 
quantities 5° and three distributions G". Keeping in mind the absence of transition of a quark 
(antiquark) to antiquark (quark) P^q = Pgg = in lowest order of perturbation theory, the evolution 
equations of these 3 sets of distributions will have a form similar to ones for quark and gluon 
densities inside a proton given above. Similar consideration was used in frames of QED in paper 
[2ll |. For quark densities 



dt ^ ' ^ 27T J y 



D%y,t)Pgg [ - ] + D^y,t)Pgg [ - 



dt ^ ' ^ 2n J y 



D''{y,t)Pgg - +DSiy,t)Pgg - 



at 271 J y 



D''iy,t)P„ 



(Al) 



Similar set with the replacement — > Z?" take place for anti-quark densities. For gluon densities 
we have 



dt 



G''iy,t)Pg, 



G\y,t)Pg, 



dt 2n J y 



G^y,t)Pgg - +G'^iy,t)P,g - +G«(y,i)Pg, - 



(A2) 



It follows from these sets of equations 



dt 



{D3{x,t) + G3{x,t)) 



a{t) f dy 
27: J y 



{Da{y,t) + Ga{y,t))P,^ 



{D^{y,t) + G^{y,t))P, 



gq 



iD''{y,t)+G^{y,t))Pg. 



dt 



{D^x,t) + G'^{x,t)) 



ajt) f dy 
271 J y 



iD''iy,t) + G\y,t))Pg, 
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+ {D<>{y,t) + G<'{y,t))P, 

When omitting the densities of anti-quarks and inside the quark and the gluon and identi- 
fying 

D'i + G" = q, 

Ds + G^ = G, 

we reproduce AltareUi-Parisi equations ([1]). 

Our statement about numerical smallness of the contribution of the intermediate quark (anti- 
quark) states in the evolution of gluon density follows from iteration procedure in solving the first 
equation of gluon set. So it can be taken into account by including as a relevant contribution to 
iC-factor. Besides only light quarks must be considered, describing the experiments without quark 
jets production. 




Appendix B: Approximate evolution equations 

A quark dominance consist in suggestion I?' ^ D'^, C. Gluon dominance imply G^ ^ D^, 
and besides ^ D^. Set of equations in these approximations reads as 

1 



X 

dt 271 J y \y 



dt 271 J y \y 



G%x,t) ^ -G^{x,t) 



dt ^ ' dt ' ' 2tt 



J y \y 



X 

Note that the equation for D'^ coincide with equation for non-singlet quark density Qns = Q ^ Q- 
The equations for D'' and G^ are given above. 

Keeping in mind the solution of evolution equations 

|a(x,/?,,) = + hpl, + ...] = ^ [a/3,, + hpl, + ...] ; 

A(x,/3,) = a^Pl+b^p^^ + ..., (B2) 



we obtain 



Da{x,l3g)^G''{x,f3,) = i(l + (l-x)2)/32 + 0(/33). 
DHx,P,)^D'''{x,(3,)^D'^'{x,f3,) - i0(x)/32 + 0(/33), 



(/,(x) = — (l-a;)(4 + 7x + 4x2) + 2(l + a;)ln2;. (B3) 
3a; 
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Appendix C: Parton densities in the proton 

Keeping in mind the problem of description of inelastic processes in high energy collision of 
protons it seems to be naturally consider proton as an objects with definite contents from quarks 
and giuons. It implied the presence of the preliminary evolution from mass shell to virtuality of 
order 1 GeV'^ of all the constituents of proton. 

Note that due to condition xiX2 > ztt only valence quark and giuons inside proton take part in 
Drell-Yan process. We will choice the density of the valence quarks and giuons approximately as 
ones found in paper fl5l |: 



xWu {x) 


= Aux^^{l-x)^-{l + euy 


+ 






= 0.2; ?7i = -0.73; m = 


3.3; 




xWd{x) 


= Adx^'{l-x)^'{l + ea^ 


fx + 




Ad 


= 18; 773 = 0, 1; 774 = 6; 






xWg{x) 


= AgX'^{l-X)^<^{l + tg^ 


/x + 


-1gX)+Ag,X^^'{l 


A, 


= 0.0012216. 







(CI) 

Numerical constants must be chosen in such a way to satisfy the constrains from number sum rules 

1 1 
dxWu{x) =2; j dxWd{x) = 1, (C2) 


and besides the momentum sum rule 

1 

dxx [W^{x) + Wd{x) + Wg[x) + S{x)] = 1, (C3) 
with S{x) is the sea contribution. 

Where 7„ = 8.9924, = 7.4730, ?7g = 2.3882, 77,,. = 0, Ag, = 0, (5g = -0.83657, 5g, = 0, 
7g — 1445.5, Eu = —2.3737, — —4.3654, eg — —38.997. More complicated expressions for 
densities which was extracted from description of fixed target HERA and Tevatron experiments 
are presented in p^ . 



R. P. Feynman, Photon-hadron interactions (Benjamin, New- York, 1972). 
H. D. Politzer, Phys. Kept. 14, 129 (1974). 

J. C. Co llins, D. E. Soper, and G. F. Sterman, Adv. Ser. Direct. High Energy Phys. 5, 1 (1988), 
|arXiv:he p-ph /04093 1 3 

C. F. von Weizsacker, Z. Phys. 88, 612 (1934). 

E. J. Williams, Phys. Rev. 45, 729 (1934). 

P. Kessler, Nuovo Cimento, X. Ser. 17, 809 (1960). 

V. N. Baier, V. S. Fadin, and V. A. Khoze, Nucl. Phys. B65, 381 (1973). 

G. AltareUi and G. Parisi, Nucl. Phys. B126, 298 (1977). 

L. N. Lipatov, Sov. J. Nucl. Phys. 20, 94 (1975), [Yad.Fiz. 20:181-198, 1974]. 

E. Kokoulina, Acta Phys. Polon. B35, 295 (2004), arXiv:hep-ph/0 401223 

E. S. Kokoulina and V. A. Nikitin, (2005), arXiv:hep-ph/0502224j 

P. F. Ermolov et al, (2005), arXiv:hep-ph/0503254 

A. Bukhvostov, G. Frolov, L. Lipatov, and E. Kuraev, Nucl.Phys. B258, 601 (1985). 

E. A. Kuraev and V. S. Fadin, Sov. J. Nucl. Phys. 41, 466 (1985). 

A. Martin, W. Stirling, R. Thorne, and G. Watt, Eur.Phys.J. C63, 189 (2009). arXiv:0901.0002] 
J. Kodaira and L. Trentadue, Phys. Lett. B112, 66 (1982). 
Particle Data Group, K. Nakamura et al, J. Phys. G37, 075021 (2010). 
N. Kidonakis, JHEP 05, Oil (2005), arXiv:hep-ph/0412422 



N. Kidonakis, Phys.Rev. D64, 014009 (2001), arXiv:hep-ph/0010002 



N. Kidonakis and R. Vogt, Phys.Rev. D68, 114014 (2003), ar Xiv:hep -ph/0308222 
A. B. Arbuzov, V. V. Bytev, E. A. Kuraev, E. Tomasi-Gustafsson, and Y. M. Bystritskiy, Phys. Part. 
Nucl. 41, 394 (2010). 

[22] F. Abe et al.. The CDF Collaboration, Phys. Rev. D 56, 2532, (1997), FERMILAB-PUB-97/093-E. 



